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Abstract 

The structure function of a heavy quark in a heavy hadron in electroproduction 
is calculated from QCD. Its deviation from the parton model prediction (a delta 
function) is shown to be governed by a shape function, similar to the one recently 
proposed by Neubert and Bigi et al. to describe the end-point behaviour of the 
lepton spectrum in semileptonic B X^ev decays and the line shape m B ^ Xg'^. 
Two new sum rules are derived which extend the Bjorken-Dunietz-Taron sum rule 
by including corrections suppressed by the inverse heavy quark mass. 



^ Supported by the GraduiertenkoUeg Teilchenphysik, Universitat Mainz 



1. There have been recently significant advances in the theoretical understanding 
of the inclusive decays of hadrons containing one heavy quark. It has been shown 
that a number of model-independent predictions can be obtained, for example for 
the lepton energy spectrum in semileptonic decays of B mesons or the photon en- 
ergy spectrum in the decay B Xg'-f |T|, |^. The formal tool which was used is 
an adaptation of the operator product expansion method (OPE) to deal with this 
specific situation in which the large mass scale (the heavy quark mass) is trans- 
ferred from the matrix elements of the operators to the Wilson coefficients. From an 
intuitive point of view, the results thus obtained can be interpreted as corrections 
to the naive free-quark decay predictions, due to a "Fermi motion" of the heavy 
quark inside the hadron. Initially introduced by Neubert to describe the end-point 
region behaviour of the lepton spectrum in B ^ XuCv [^], the concept of a "shape 
function" turned out to be a process-independent parametrization of the "Fermi mo- 
tion" , incorporating features which cannot be reproduced by a quantum-mechanical 
description in terms of a wavefunction @, |^, 0] • 

In this paper we discuss the h quark structure functions in electroproduction on a 
B meson and show that they can be calculated in QCD into an expansion in powers 
of A/mfe. We give explicit results for these function to second order in this parameter. 
To leading order in this expansion the parton model prediction is obtained, which is 
a simple delta function. The corrections to it have a singular form consisting of delta 
functions and its derivatives. It is shown that the most singular contributions can 
be summed up to all orders in l/m;, into a "shape function" which extends beyond 
the parton model region and which can be considered a good approximation to the 
real distribution function. 

In Section 4, two generalizations of the Bjorken sum rule are derived [Q. The 
sum rules result from equating the total contribution to the moments of the heavy 
quark structure functions in electroproduction on a heavy meson (as obtained by 
the OPE analysis) to the sum of the individual contributions of the resonances. 
Both these quantities can be arranged as an expansion in l/mq, the inverse heavy 
quark mass, and at the leading order in this expansion the usual Bjorken sum rule 
is recovered. At higher orders in l/mg the resulting sum rules are new and can be 
used to set constraints on the nonperturbative subleading form-factors. 

The spirit of the approach and the nature of the approximations which are made 
are related to those recently used [0 to calculate the inclusive semileptonic decays 
of heavy hadrons containing a heavy quark. However, the formalism used here is 
more similar to the one familiar from deep-inelastic scattering. Instead of modifying 
the operator-product expansion (OPE) so that the heavy-mass dependence resides 
from beginning in the coefficient functions, we will prefer to perform the OPE in 
the full theory (QCD), and to expand the resulting matrix elements in powers of 
l/mq. The equivalence of the two approaches is proved on an explicit example to 
order l/mq. 
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The heavy quark structure functions in electroproduction on a heavy meson have 
been considered previously by Jaffe and Randall 0. In this paper we go further by 
giving explicit expressions for quantities which were just parametrized in 0. Also, 
some of the results described in Section 3 have been obtained by Neubert in a 
light-cone language. Our derivation is perhaps more direct as it proceeds directly 
from the usual definitions of the structure functions. 



2. Let us consider the forward Compton scattering matrix element on a heavy 
meson B 

T,,iq-P,Q') = -ij ^'xe-''^-%B{v)\T{hM^){hMmB{v)) , (1) 

where p = itlbV is the B-meson momentum and = —q^ > is the photon vir- 
tuality. The meson state \B{y)) has the usual normalization 2mBVo- This situation 
is similar to what one usually has in deep-inelastic scattering, in that the virtual 
photon is spacelike. However, in distinction to the latter, we will not require that 

^ m^, but will take Q^, m\ ^ ^qcd (with the exception of Section 3, where we 
take temporarily oo). We will keep therefore all contribution of order m?^/Q'^ 

and organize the result as an expansion in powers of A/m^. 

T^y can be decomposed into covariants defined in the usual way by 

T,M-p.Q^)= (2) 

-Ti{g,,u - + ^(p^ - q^.^){Pu - qu-^) - ^T^ie ^.^pp'' q" . 

In our case, because the electromagnetic interactions are parity-invariant, T3 = 0. 
The analytic structure of the scalar form-factors Ti ■ p, Q"^) in the q ■ p complex 
plane for fixed is shown in Fig.l. It has two poles ai q-p = ±(5^/2, corresponding 
to the elastic scattering process 7*5 ■y*B and two cuts for \q-p\ > l/2[Q'^ + {mB + 
m^)^ — m^]. Ti^2 are even functions in the variable q ■ p: 

Ti,2(-g-P,Q') = Ti,2(g-p,Q'). (3) 

The discontinuities of Ti^2 across the cuts are related to the hadronic tensor Wfj,^ 
which describes the inclusive process B ^ e^Xh (in the approximation that only 
the b-quark is struck): 

= ^T.iB\ihMm){X\{hubm\B)i2nr6ipx-p + q) (4) 
47r ^ 

which has a decomposition into covariants similar to T^^. We have 

discT^j,(g,p) = AmW^.^ forf-g<0 (5) 
discT^i.(g,p) = -discT^y{-q,p) fort;-g>0. (6) 
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To identify the operators which appear in the OPE (1), the usual method is to 
take the matrix element of the time-ordered product between free 6-quark states. 
This gives 



d^xe-"^-^b{r)\T{b^^b){x){b^,b){0)\b{r)) = (7) 

^ i/-</+mb f ^ , -f ^ !/+</+ nib , X 

"0" 7^7 2^Mr) + u{r)-f^ 2^M'^) 

[r — qY — mi [r + qY — mf, 

with r the momentum of the quark. This has to be expanded in powers of r and 
the result identified with the matrix elements of a sequence of operators of the type 
b - ■ - b with an increasing number of derivatives acting on the fermion field (one 
for each power of r). At leading order in the coupling (7, we can replace in the 
denominators of (7) by m^, since it generates contributions of the form {iDYb which 
equals (in matrix elements) mfb in virtue of the equation of motion for the b field. 
The final result is 

-t J d^xe-*'?-^T(67^6)(x)(67.&)(0) = (8) 

~Q4 ^ [9fii-ii9u^i2 ~ 1^9 nu9 ^J.l^l2\ ~^ QfiQfii9i'fi2 ~^ Q^ilQl^9^l^l2 ~ QtJiiQiJ.29 iJ.u} 

2m °° / 2 \ 

n ■■■n pfJ.llJ.2-- -1^271+2 tn \^ I 1 n ■■■n UfJilJ.2-fJ.2n 

Hfl3 HfJ.2n + 2 Q2 ilfJf X^Q"^ j ^fJ.2n-^^ 

Here P and R are defined to be the operators 

Pfj.ifJ2-fin = ^ (&7mi(^-D)m2 ■■■ (^-D)Mn^ + permutations of the indices) (9) 
QfJifJ2-fJn = ^ •■• + permutations of the indices) . (10) 

In deriving (8) all operators but the completely symmetric ones have been neglected, 
as their matrix elements in an unpolarized hadronic state vanish. This is the relation 
on which all our discussion will be based. We will take its expectation value between 
heavy hadron states and the matrix elements of the operators P and R will be 
expressed as power series in 1/mf, by using heavy quark effective theory methods. 
But before doing so, we pause for a moment to discuss the method by which the 
information obtained with the help of the OPE is related to physical quantities. 

Previous studies of the inclusive semileptonic decays of hadrons containing 
one heavy quark have simply identified the hadronic tensor with the discontinuity of 
the amplitude T^j^ calculated using perturbative QCD. The results thus obtained are 
highly singular and must be understood in the sense of duality, when integrated with 
an arbitrary smooth function. In this paper we will take a different approach and 
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will relate the moments of the structure function to the coefficients ofv-q appearing 
into an expansion of T^^ around the origin in the v ■ g-complex plane with the help 



of a dispersion relation |Ty]. In the end the results obtained with the two methods 
will be seen to be identical. The principle of the method is the following: for very 
large Q^, the two cuts on Fig.l are well separated and one can use perturbative 
QCD to reliably compute T^i, around the point q ■ p=0. In this region each of the 
invariant functions Ti^2 can be expanded in a power series 

T.iq-p,Q') = f:TUQ')(^^J ■ (11) 
The coefficients can be expressed as a contour integral 

T^niQ') = ^ I ck%^ (12) 



2T[i Jc z""*"^ 

where the contour C is a small circle around the origin. It can be deformed to the 
shape shown in Fig.l, where it picks contributions from the discontinuities across 
the cuts. The contribution of the large circle at infinity will be neglected. By making 
use of (5,6), the contour integral in (12) can be written as 

TiniQ^) = -2 /' d>:x^~'W,{x, Q2) (13) 
Jo 

where x = Q^/(2g ■ p) is the Bjorken variable. This is the desired relation which 
expresses the moments of the hadronic tensor in terms of quantities calculable the- 
oretically, the expansion coefficients Tj„. 

We turn now to the evaluation of the transition matrix element T^^, by taking 
the matrix element of the OPE expansion (8) between heavy meson states. This 
leads us to the problem of calculating the matrix elements 

\B{v)) {B{v)\R,, \B{v)) . (14) 

They can be expressed as power series in the inverse heavy quark mass 1/mf, by 
making use of heavy quark effective theory methods. First, the b-quark field is 
written as 

b{x) = e-'"''''-'=B{x) (15) 

where 



ip^ 1 / 1 



The heavy quark field h^^ satisfies i/h^^ = h!;^\ 



h^^\x). (16) 
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Let us consider for beginning the matrix element of P. For simplicity we exem- 
plify the calculation on the first term in (9), the final result having to be obtained 
by symmetrization with respect to all indices. First we write 

{B(v)\h,AiD),, ■ ■ ■ (iD),J\B{v)) = (17) 
{B{v)\B-fu_^{mbV + iD)^^{mi,v + iD)^^ ■ ■ ■ {rubV + iD)^^B\B{v)) . 

To leading order in the large mass scale this is equal to 

2mB mp^-y^i-y^j • • • "^^in , (18) 

which is already symmetric in all indices. Analogously, the matrix element of R can 
be written as 

{B{v)\B{rrn,v + iD)^,{mi,v + iD)^, ■ ■ ■ {m^,v + %D)^^B\B{v)) = (19) 
2mBm^Vn^Vn^ ■ ■ ■ v^^ + terms suppressed by l/mi,. 

If we insert these expressions for the matrix elements in the equation (8) for T^,^, 
the following results are obtained for the invariant transition amplitudes 



Ti{q-p,Q') = - 
T2{q-P,Q') = - 




If we define the more conventional structure functions Fi{x,Q^) = Wi{x,Q^) and 
F2{x,Q'^) = q • plm\W2{x,Q'^), the following predictions are obtained for their 
moments, with the help of Eq.(13) 

/^dra;'^-^Fi(a;,g2) 

^0 

/'drx"-^F2(a;,g2) 

^0 

The evenness property of 71,2 (3) prevents us from extracting the moments (22,23) 
for other vahics of n than those indicated. However, these can be obtained by 
analytic continuation, that is, by simply assuming that these relations hold true for 
all values of n. 

Some comments can be made about this result: i) the integral of Fi(a;, Q^) 
equals unity (n=l), which corresponds to the parton interpretation of this function 
as being the probability for the heavy quark to be found with a fraction x of the 



n— 1 



/mb_ 
\mB 

KitibJ 



1 - 



2 1- 



A_^ 

TUB 



rriB, 



n-l 



n = 2,4,6,... (22) 
n = 1,3,5,... (23) 
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total momentum; ii) as expected, the Callan-Gross relation Fl = F2 — 2xFi = is 
satisfied; iii) the shape of the distribution function is a simple delta-function 



F,ix,Q') = 6ix-^) (24) 
rriB 

which is the naive parton model result, although modified to take into account the 
binding energy of the heavy quark in the hadron. 

Let us now go further and examine the effect of retaining more terms in the 
matrix elements (17) and (19), suppressed by one or more powers of l/ruf,. The 
correction of order l/nib to (17) is equal to 



ruf, 



n-l 



E^M2 ■ ■ ■ ■ ■ ■v,AB{vM\AtD),M';^^\B{v)) (25) 



i=2 



n -if) 
j=2 ^ 

j=2 ^ 

■j ^ r - 



j=2 

The matrix element in the first term can be written as 

{Biv)\h^:h,.m,X:^\Biv)) = Av,,v,^ (26) 

with 

A = {B{v)\h'^\iv ■ D)h^^^\B{v)) = (27) 

where the equation of motion iv ■ Dh^^ = has been used. The other three terms 
have together the form of the l/m^ correction to the matrix element of the vector 
current &7^i&, which is known to have no such corrections. Therefore the entire l/mb 
correction to the matrix element of the P operator vanishes. A similar argument 
shows that the l/m^ corrections to the matrix elements of R vanish also. From 
this follows the conclusion that there are no l/rub corrections to the distribution 
functions of a heavy quark in a heavy hadron. 

In principle the same line of reasoning will give also the 1/m^ corrections to the 
distribution functions. This will require the introduction in the operator product 
expansion (8) of some new twist three operators containing one gluon field tensor 
F^i/. Much simpler it is however to use the result for T^^ given in the second paper 
of Ref. nil. We will show that the answer obtained with the two methods for the 
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twist-2 contribution is the same. In the next section we discuss the corrections of 
order 1 / ml to the distribution functions of the heavy quarks. 



3. Blok, Koyrakh, Shifman and Vainshtein have calculated T^^, up to order 
A'^/ml (Eqs.(Al,A2) in the second paper of Ref. Although they were interested 
in a timelike q, their expressions are also valid for a spacelike q, the case we need 
here. The results for Ti^2 are 

^ .T, = h-vq + ^{K, + G,)] (28) 



2mB A' 3 
2m> 1 

+ ^{-^Gk[-iv ■ g)(4m, - 3v q) + 2[{v ■ qf - q^ 
+K,[-{vqf-\[{v.qf-q']]} 

^'^'K,{vq)[{vqf-q'] + {q^-q) 



3A3 



and 



1 Q,Tn 1, ^ Ath'^ 

-T2 = -r^[l + 7t(^6 + G,)] + -ri(7ir, + 5G,)v ■ q 



with 



2mB A ' 3' " 3A^ 

+ ^-^K,[{v-qr-q'] + {q^-q) (29) 



A = -2mbV -q-Q^ + ie. (30) 



The terms with q —q have to be added to account for the necessary crossing 
property (3) of Ti 2- Kf, and Gh are defined by 



. _J_(B,.„)|/4'.(||.f|B,.„)) ^ ^ ^ (31) 



and have the approximative values Gb = —0.0077 and Kb = 0.01. We have shown 
here the different names by which these quantities are denoted in the literature 

B i MM- 

The expressions (28,29) are the first terms of a power series in the expansion 
parameters A/rrib and A^/A. Because we are interested in evaluating these functions 
around the point v ■ q = 0, the second expansion parameter is actually K^/Q^. 
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We expand the expressions for Ti^2 into a power series in {2mBV • qj/Q"^ with the 
result for the coefficient of the n*^ power 

T,„.-2(^)"-'{l + (n-l,5^^±i^ 

and 



(34) 



for n even (0,2,4,...) and Tj „ = for n odd. 

According to the discussion in the preceding Section, these coefficients are di- 
rectly related to the (n — Vf^ moment of the respective structure functions. By 
performing an inverse Mellin transform or simply by trial, the functional depen- 
dence of F\^2{x).i including these corrections, turns out to be of the form 

V 3 Q'^ J TUB 

nib f 5 ^^ c^n^ 2 m?\ , mb 



-{Kb + Gb) + S{Kb + -G,)-| 5'{x - 
ruB \o 6 Q'' J ruB 



rrib 



3 KmB/ \ I ruB 



and 



\ (^) ' xF,{x, Q^) = f 1 + AKb + ^G, + SKb"^) 5{x - ^) (36) 



3 Vms/ \ Q^/ niB 

This is the same as what one obtains by directly taking the discontinuity of 
(28,29) by using the well-known identity 
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The Eqs. (35,36) are all what QCD predicts about the structure functions of a 
heavy quark in a B meson at order A^/m^,A^/Q^. They have a rather singular 
dependence of the Bjorken variable x, given by a 5 function and its derivatives. 
These singularities are similar to what has been recently obtained for the shape of 
the electron spectrum near its end-point in i? X^eu or for the photon spectrum 
in B ^ Xsj decays. An interesting feature of these predictions is that, although the 
support of the delta functions is restrained to a single point x = mb/mB, they do 
simulate a spreading of the structure functions from the parton-model prediction. 
For example, one can compute the dispersion of the x variable for the structure 
function Fi{x): 

al = {x'),, - {x)l^ = ^ (^IXb - S{Kb + ^G,)^) (38) 

which is nonvanishing and is of the order A/mt. 

Of course, the physical structure function does not look like a delta function and 
derivatives thereof, but has a smooth shape which extends beyond the parton model 
boundary^. We will show that the sum of the most singular terms in the expansions 
(35,36) to all orders in A/mt generates a "shape function" which can be considered 
as an approximation to the physical structure function. The argument goes in close 
analogy to the analysis of Neubert of the photon spectrum in B ^ Xg'y decays. 
However, the details of the derivation are somewhat different, as appropriate to the 
problem at hand. 

We will consider in the following the scaling limit oo when the distri- 

bution functions depend only on the Bjorken variable x. Note that these are not 
the distribution functions which would be measured in a hypothetical electropro- 
duction experiment on B mesons, since the relevant matrix elements are evaluated 
at the heavy quark scale and the observable distribution functions would have to 
be obtained by evolution to the large momentum scale Q^. However, we will have 
nothing to say about this evolution and therefore all our results refer exclusively to 
the "boundary data" for the distribution functions at the scale /i^ = m^. 

In this limit there is only one independent structure function, since the two 
functions Fi{x) and ^2(0:) are connected by the Callan-Gross relation F2 = 2xFi. 
One can easily convince himself that the predictions (35,36) satisfy this equality in 
the limit Q'^ —>■ 00. It is only violated at first order in as- Therefore, we will proceed 
further with Fi (x) and call it simply "the distribution function" . 

The physical distribution function in the scaling limit Fi{x) can be represented, 
when convoluted with an arbitrary smooth function, as the sum of a delta function 

^We assume that the rapid variations in the physical structure functions due to the resonances 
have been smoothed out. 



9 



and its derivatives: 



i^i(^)=f:(-ir^<^«(^-^) (39) 



where the coefficients Mj are the moments of Fi[x) 

Mi = C - —YFAx) . (40) 

JO TJlB 

The first few moments Mj can be extracted from Eq.(35) 

Mo = 1 (41) 
= ^^(i^, + G,) (42) 

6 TUB 

M2 = l^Kb (43) 
3 

Each of the moments Mj can be expanded in a power series in A/m^ and it is 
reasonable to assume that keeping only the leading term in this expansion provides 
a good first approximation to the description of Fi{x). This is equivalent to keeping 
only the leading singularities in (38) (the smallest power of A/nib corresponding to 
a derivative of a given order of the S function is the same as the highest order of the 
derivative corresponding to a given order in A/mi). 

In the scaling limit ^ oo only the operator P in (9) gives a nonvanishing 
contribution to the distribution functions, since it is the only one which contains 
twist-2 operators. The operator it! and any other operators involving the gluon field 
tensor give contributions suppressed by one or more powers of Q^. The matrix 
elements of P have the form 

iM2 -M.J-^('^)) ~ C'^tii'^ti2 ■ ■ ■'^/in + terms containing g^yS (44) 

We are only interested in the part proportional to v^s because the terms with g^^ 
vanish in the scaling limit. Let us try to extract as much information as possible 
about the constant C to all orders in l/rrib- It has an expansion of the form 

C = IrriBml-^ (l + eCi + e^C^ + • • •) (45) 

with 

6 = — (46) 
mb 

A direct calculation by making use of the methods exposed in Section 2 gives 

eCi = (47) 
e'C2 = ^X,^^^:^^^^^ + ^(X, + G',)(n-1) + (X, + G,)(n-1) (48) 
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The first relation expresses tlie vanisliing of tlie l/m^ corrections to tlie structure 
functions, noted in the preceding Section. The three terms in the second relation 
have respectively, the following origin: i) the first appears from the matrix element 
{B {v)\{iD) fj^{iD) {v)) resulting from the expansion of (17); ii) the second appears 
from the correction to the B field (16) at order l/mf, and iii) the third appears from 
the insertion of the l/rUb term in the HQET Lagrangian. 

If the expression (44) together with (45,47,48) are introduced into Eq.(8) and 
the result for Ti expanded in a power series in {2mBV • q)/Q'^.i the coefficient of the 
n*'* power which is obtained coincides with Eq.(33) (in the hmit oo). This 

establishes the equivalence of the two approaches to this order in A/m;,. 

The constant C is directly related to the (n — 1)*'' moment of the structure 
function Fi(x), as can be seen upon inserting (45) into the general formula (8): 



L 



drx"-^Fi (x) = V" ^ (l + eCi + e'C2 + • • •) (49) 



10 XmB/ 

For a given order in e (say, e'^), the most singular terms in Fi correspond to those 
terms in the coefficient Ck which have the highest power of n. For example, the term 
with the highest power of n in C2 is ^Khu'^ and it corresponds to the contribution 

(^) ^"(^ ~ ^) -^i(^)) Eq.(35). As said before, this term originated from 
simply expanding the brackets in 

{B{v)\B{mbV + iD)^, {mbV + iD)^, ■ ■ ■ {m^v + iD)^„B\B{v)) , (50) 

replacing B by the HQET field h^'^ and forgetting about lagrangian insertions. We 
will show now by a simple counting argument that this holds true generally and that 
the leading terms in n are correctly reproduced by retaining only these contributions, 
to all orders in e. 

The contribution to of these terms can be written as 

Y.{B{v)\h^h,AiDUme.. ■ ■ ■ (iD)^M^\B{v))v^^^, ---v^^ (51) 

where the sum runs over all possible ways to choose the k derivatives out of a total 
of (n — 1) possibilities, and contains thus terms. 

Let us consider now a contribution to Ck arising from corrections to the B field. 
We write the expansion (16) of S in powers of l/m;, in a compact form as 

00 ry 
n=0 

with Bq — l^Bi — {ip)/2, etc. Then the desired contribution reads 

5: {B{vM^Ba,,m^,{iDU ■ ■ ■ m^^_^_^B,h'^^\B{v))vo.,_,_^^, 
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and the sum contains 

k-lk-i-l 1^^^ _ ^y^ fc-1 ^. ^ 

terms. We can omit the terms with i + j = k because this is a correction to the 
matrix element of the vector current which is known to vanish. We have included 
here also the preceding case as the term with i = j = respectively r = i + j = 0. 
Assuming all the matrix elements in (53) to have about the same order of magnitude, 
one can see that the total contribution from the matrix elements without corrections 
to the B field (r = 0) scales hke n'' for a given order k in the e expansion. The 
next contributions, containing a suppression by l/m^ from a correction to the heavy 
quark field (r = 1), scale as n^~^, and so on. The matrix elements containing 
insertions of higher-dimensional terms in the HQET lagrangian can be shown in a 
similar fashion to give only subleading contributions in n. This proves the above- 
mentioned assertion that the most singular terms in Fi{x) are correctly reproduced 
by replacing 

{B{v)\P,,,,...,jB{v)) ^ (55) 
Vf,,{B{v)\U^\mbV + iD)^,{mbV + iD)^, ■ ■ ■ {m^v + iD)^^h^^^\B{v)) 
+ symmetrization w.r.t. all indices 

We define the dimensionless constants a„ ~ 1 through 
1 



B{v)\h^\iD),,{iD),, ■ ■ . {iD),^h^)\B{v)) = a^AX.^M. • • "^^^ + • • • ■ (56) 



2m B 

The first few Oj are 

ao = 1 (57) 
eai = (58) 

= ^Xfe. (59) 

In terms of these constants the matrix element (55) can be written as 
v^,{B{v)\U^\mbV + iD)^,{mbV + iD)^, ■ ■ ■ {rrn,v + iD)^^h^^^\B{v)) (60) 

.1-1. . „-i^' (n- 1)! 

+ symmetrization w.r.t. all indices = zmBirij^ 2^ tt-, : rraj — jV^^Vn^ ■ ■ ■ v^^ . 

i=o i'-[n — i — 1)1 ml 

Inserting this in Eq.(8) gives for the moments of Fi{x) (only the part containing the 
leading singularities, the so-called "shape function" ) 
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The moments Mi with respect to the parton- model point x = rrib/mB which 
appear in (39) can be obtained from the above ones by a simple expansion 



Jo niB ^«!(n-z-l)!V mB/ Jo 

= r^V'^ V vr-n— 1 = r!!^^ 



(62) 

where the final result has been obtained by interchanging the summation order. 
It has been shown |^ that the line shape in the B Xgj decays and the 

end-point electron spectrum in the semileptonic B — * Xu decays are essentially 
determined by one common universal function whose moments are related to the 
nonperturbative matrix elements in a fashion similar to (62). Therefore the result 
(62) is remarkable, as it shows that the same universal function gives also the heavy 
quark structure functions for electroproduction on a B meson. 

Most probably, these structure functions will never be measured and therefore 
our results are of little practical value. Nevertheless, we believe that they are in- 
teresting in themselves because they show how completely different phenomena in- 
volving heavy quarks are governed by similar distribution functions. Finally, we 
mention the related problem, of considerable experimental interest, of the fragmen- 
tation functions for heavy quarks into heavy hadrons, for which one cannot make 
use of the powerful formalism of the operator product expansion. One can hope 
that similar methods could be used to study the shape of these functions, and set 
constraints on the possible structure of their nonperturbative contributions. 



4. In this Section we derive two new sum rules for the nonperturbative functions 
which describe the form-factors for heavy-to-heavy transitions. These sum rules are 
deduced for the process •y* + B ^ Xf,, with Xf, an arbitrary excited state of the B 
meson, but they can be extended by using heavy quark symmetry for the case of 
greater interest B ^ Xc + W, with Xc an excited state of the D meson. 

We will relax the assumption — > oo and keep arbitrary (but still much 
larger than Aqcd)- Let us examine the form of the physical structure functions 
Fi 2(a;, Q^). In the neighbourhood of x = 1 (near the end of the cuts) they have a 
rapid variation, with spikes corresponding to each resonance Xf, superposed over a 
smooth background due to two- and many-body intermediate states, such as Bn. 
The contribution of a particular resonance of mass mx to W^u can be obtained from 
(4) to be given by 

Wl,^\p ■ q, Q') = ^6[p . g - l(m| - - Q^)] x (63) 
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{B{v) I {h,h) (0) \X,{v')){X,{v') I (67.6) mB{v)) 
where the velocity v' is given by 

f ■ w = . (64) 

The matrix elements can be expanded in powers of l/m^ by using heavy quark 
effective theory methods and the resulting structure functions Fi 2 are expressed in 
terms of the usual Bjorken variable x by using 

5[p-q- -{ml -m\- Q^)\ = -^6{x - xx) (65) 



with 

The key ingredient which makes possible writing a sum rule is the positivity 



Xx = -p^ :7i— — 2~ • (66) 



property of F12 ||T2[) which expresses the fact that each individual intermediate 
state in (4) contributes a positive quantity. Equating the OPE result for the mo- 
ments (33,34) with the moments of the structure functions written as sums over 
intermediate states gives an upper bound on any partial sum of the latter. 

On the tip of the cut, mx = ^ub and v-v' = 1 + {2m1). The velocity transfer 
V ■ v' decreases as we move along the cut until it reaches a minimum [v ■ v')min = 



1 + Q"^ /m\ for nix = f^By^ + Q"^ I'^^'b after which it starts increasing. Therefore 
the various intermediate states will appear in the sum rule with a different velocity 
transfer v ■ v' . 

From (63) we get 



^ 2[Q^ + rrix — mgY 

(i?(^)|(M)(0)|X,(t;'))(^.(^')l(M)(0)|i?(t;)) 

or, in explicit form for the s- and p-wave mesons 

Q') = nfn2^ T ^^(^ - ^xjR' (68) 

with 

• B meson R = 0. 

• B* meson R = niBmB^iu^ - 1)1^ + i;(2xi - 2(cu - l)x2 + 4x3 + - 6)1'- 
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• Bo meson (s^^ 


= 1/2+ , s-^ 


= 0+) R 


= 0. 






• Bi meson {s^ 


= l/2+,s^ 


= 1+) R 


= \mBmB^ 


{u- 


- mil/2?. 


• B[ meson (sj* 


= 3/2+, 


= 1+) R 


= ImBiriB'^ 






• B2 meson (sj* 


= 3/2+, 


= 2+) R 


= \mBmB^ 






and for xF2{x, Q^) 













(g2 + mV - m|)2 



(69) 



with 



B meson S" = |^ + ^(xi - 2{uj - l)x2 + 6x3 



B* meson S = 

M-is)?. 

Bq meson (sj* 
Bi meson (s^' 



4mB 



4 - 



4 



+ ) le + i;(2xi - - 1)X2 + 4X3 + 



1/2+, = 0+) S = 0. 



1/2 



B[ meson (s 

1)|6/2P. 

B2 meson (sV 



I 



3/2- 



1+) 5 =1(^-1)1^1/2 

= 1+) S 



8mB 



1 + ^ 



+ 



2mB 
3m „/ 



(^ + 4) (^2 _ 



3/2+, s'^ = 2+) S 



SniB 



-B2 



+ 2^-)(-'-l)le3/2p. 



Here u stands for the velocity change v ■ v' . ^(tu), ^i/2,3/2(^^) are the Isgur-Wise 
functions which describe the form-factors of the transitions from the respective mul- 
tiplets to the lowest-lying one0. For the s-wave multiplet we have included also 
the contribution from the subleading form-factors Xi, 2,3,^8! which are defined as in 



Tj, ITJ]. For example, if we assume for the moment that the multiplets are degen- 



erate, their contribution to Fi^2 takes the form 
Fr(x,Q2) 



+ —(2X1 - 2{u - 1)X2 + 4X3 + M- i3)?S{x - 1) 
2 nil, 



Q^[Q' + {mB,-mBf 

6[Q2 + «^-m|)p 

Q^[Q^ + [mB,-mBf 
3[g2 + (^|^_^| 



2 M2 



■{u - l)\iii2?5{x - XBi) 

iuJ-l){l + Uj)\3/2?6{x-XB,) + 



(70) 



^^1/2,3/2 are related to the form-factors used in Q by ^1/2 = 2\/3ti/2 and ^3/2 = V3t3/2- 
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and 



xFr{x,Q')=^-^\iMx-l) 



Q 



4 



6[g2 + (m|^ - ml)f 



u:-l)\iy,YKx-XB,) (71) 



Q 



4 



3[g2 + (^1^ _ ^|)]2 



(u;-l)(l + a;)^|e3/2r(5(x-XBj + 



However, since we work to an arbitrary order in l/rrib, wc will consider henceforth 
the masses of the two members of each multiplet as different. 

The sum rules result when the moments of Fi^2 expressed as a sum over reso- 
nances are equated with the QCD prediction 



^0 KijibJ 3 

4 /m 
3 \Q 



A { ^1 + 1)^6 + {n^ + 3n + A)K,] + • • •} (72) 



and 



ymsJ 3 

+^X,(n+l)(n + 2)(^^)+---} (73) 

At this point a serious problem emerges which might diminish the usefulness of 
the sum rules. The most interesting region where one would like to apply these 
sum rules is in the vicinity of the equal- velocity point v ■ v' = 1, which corresponds 
to = 0. But, as one can see from (72,73), the l/ml corrections become very 
large in this limit, and therefore one can trust no longer the expansion in powers 
of l/mft. This happens because the expansion parameter in this problem is A/Q. 
For example, if we require as a limit of applicability that A/Q — 0.5, this would 
correspond to f ■ f ' ~ 1.018 on the tip of the cut. 

There exists yet another way to apply the sum rules, which allows one to consider 
practically any value of v ■ v'. This can be done if the b quark mass is considered as a 
free parameter. Then one can proceed in the following way: i) the ratio Q"^ /m\ must 
be kept fixed to a value given by the velocity change v ■ v' one wishes to investigate, 
according to (64); ii) and Q must be sent simultaneously to infinity (compared 
to A, which is of course fixed); iii) the result is expanded in inverse powers of m^. 
This procedure makes sure that the higher order corrections in (72,73) stay smaller 
than the leading term. 
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In the limiting case when Q'^^ms — > c>o with Q'^/m?^ = 2{uj — 1) the two sum 
rules for Fi and XF2 reduce to the usual Bjorken sum rule [ic 



1 = ^ler + - i)i6/2r + ^(.^ - i)(i + ujf\i,„\' +.... (74) 

For each n, a different sum rule is obtained. Higher values of n suppress the 
contributions of the excited states, since xx is smaller for these. We note that one 
cannot write one sum rule for each power of since the various contributions 

will not contribute in general with the same sign, so we cannot set constraints on 
the subleading form-factors appearing to a given order in the heavy mass expansion. 
Rather, the sum rules can be used to test the consistency of a set of calculated form- 
factors: the total contribution of a number of resonant states, calculated to a given 
order in l/rrift, must not exceed the QCD prediction for the sum of all excited states, 
calculated to the same order in l/m^. 

Note added. After completing this work, the paper |]T3| appeared, which extends 
the results in 0, |^, ^ ^ to cases when the final quark mass is nonvanishing. Also, 



a different derivation of the Bjorken sum rule has been recently given in |16 . 
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Figure Caption 

Fig.l The analytical structure of the forward Compton scattering amplitude T^^ as 
a function of p ■ g for given Q^. 
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